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Abstract — We investigate the stability of yr-phase in a 
polarized superfluid Fermi gas (N^ > N^, where No- is 
the number of atoms in the hyperflne state described by 
pseudospin-a). In our previous paper [T. Kashimura, S. 
Tsuchiya, and Y. Ohashi, Phys. Rev. A 82, 033617 (2010)], we 
showed that excess atoms (A A*' = Nf — N_i) localized around a 
potential barrier embedded in the system induces the vr-phase 
at T = 0, where the phase of superfluid order parameter 
differ by tt across the junction. In this paper, we extend 
our previous work to include temperature effects within the 
mean-held theory. We show that the vr-phase is stable even 
at finite temperatures, although transition from the vr-phase 
to 0-phase eventually occurs at a certain temperature. Our 
results indicate that the vr-phase is experimentally accessible 
in cold Fermi gases. 

Index Terms — cold Fermi gas, Fermi superfluid, ferromag- 
net, TT-phase, magnetic junction 



I. Introduction 

THE recently realized superfluid ^°K [1] and ^Li [2], 
[3], [4] Fermi gases are expected as useful 'quantum 
simulators' for metallic superconductivity. This is because, 
while the background physics of superfluid Fermi gases is 
the same as superconductivity, the former system has the 
unique property that various physical parameters can be 
experimentally tuned. In particular. We can study superfluid 
properties of this system from the weak-coupling regime 
to the strong-coupling regime in a unified manner [5] 
(which is sometimes referred to the BCS-BEC crossover 
in the literature), by adjusting a tunable pairing interaction 
associated with a Feshbach resonance [6]. Thus, this system 
is expected to be useful for the study of strongly correlated 
high-Tc cuprates. Indeed, using this advantage, the so- 
called pseudogap phenomenon, which has been extensively 
discussed in the underdoped regime of high-Tc cuprates 
[7], [8], was recently observed in the BCS-BEC crossover 
regime of "^^K Fermi gases [9], [10]. 

In our previous paper [11], as another application 
of superfluid Fermi gases to condensed matter physics, 
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we theoretically proposed an idea to simulate a super- 
fluid/ferromagnet/superfluid (SFS) -junction. By numeri- 
cally solving the mean-field Bogoliubov-de Gennes theory 
in real space at T = 0, we showed that, when we put 
a nonmagnetic potential barrier in a polarized superfluid 
Fermi gas (N^ > iV^, where No- is the number of 
Fermi atoms in the atomic hyperflne state described by 
pseudospin-cr), it is magnetized in the sense that some of 
excess t-spin atoms are localized around it. The superfluid 
Fermi gas is then separated by this pseudo-ferromagnetic 
junction, the structure of which is similar to a supercon- 
ductor/ferromagnet/superconductor (SC/F/SC)-junction in 
metallic superconductors [12]. We also showed that the 
SFS-junction really works as a ferromagnetic junction, be- 
cause the so-called 7r-phase (which is a typical phenomenon 
realized in the SC/F/SC-junction) is stably realized, where 
the superfluid order parameter changes its sign across the 
junction. 

In this paper, we extend our previous work at T = 
[11] to the case of finite temperatures. Since experiments 
are always done at finite temperatures, this extension is 
important to clarify to what extent the SFS-junction, as 
well as the 7r-phase, are stable against thermal effects. 
In a two-dimensional attractive Fermi Hubbard model 
with population imbalance, we examine the superfluid 
order parameter, as well as local population imbalance (or 
magnetization), by numerically solving the Bogoliubov-de 
Gennes equations. We show that the 7r-phase stably exists 
even at finite temperatures, so that one can experimentally 
access this interesting phase. 

II. Formalism 

We consider a two-component Fermi gas decribed by 
pseudospin a =t, i- To investigate the stability of vr-phase 
induced by a SFS-junction in a simple manner, we treat 
a polarized Fermi gas described by the two-dimensional 
Hubbard Hamiltonian, 
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Fig. 1 . Spatial variation of model double- well potential Vi used in this 
paper. We take V*^ = O.OOlt, V^^ = O.lt, V^ = O.lt, and ^ = 5 
(where the lattice constant is taken to be unity). These parameters are 
used throughout this paper 



Here, cj ^ is the creation operator of a Fermi atom with 
pseudospin a =t, I at the i-ih site. — t is a nearest- 
neighbour hopping energy, and the summation {i,j) is 
taken over the nearest-neighbour pairs. — [/(< 0) is an on- 
site pairing interaction, and n^cr = c] ^Ci^cr is the number 
operator at the i-th site. Since we consider a polarized 
system in this paper, the chemical potential fi^ depends 
on pseudospin a. Vi is a double-well potential, given by 



V = Vlil + V^ytl + Vo' exp [-{t./if 



(2) 



where i — {ix,iy) is the spatial position of the i-th lattice 
site. The first two terms describe a cigar trap (or anitotropic 
harmonic trap), and the last term gives a potential barrier 
around i^ — 0- i is the width of the barrier potential. We 
explicitly show the spatial variation of Eq. (2) in Fig. 1. 

In the mean-field approximation, the model Hamiltonian 
in Eq. (1) reduces to 
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where A^ = U{ci^iCi,^) is the superfluid order parameter, 
which is taken to be real in this paper. Noting that Eq. (3) is 
a bilinear form with respect to the creation and annihilation 
operators, as usual, one can conveniently diagonalize it by 
the Bogoliubov transformation [11], 

/ ci.t \ / ai,t \ 
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Here, W is a 2M x 2Af -orthogonal matrix, where M is the 
total number of lattice sites. We briefly note that the deter- 
mination of W corresponds to solving the Bogoliubov-de 
Gennes equations in the case of uniform Fermi superfluid 
[13]. 

The diagonaUzed mean-field Hamiltonian has the form 
[11] 
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where Ej^^ is the Bogoliubov single-particle excitation 
spectrum, and 
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The superfluid order parameter A^ and the particle 
density {"fii.a) are given by, respectively. 
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Here, f{E) = l/{e^/'^ + 1) is the Fermi distribution 
function (where we take fee = 1). At T = 0, the 
Fermi distribution function reduces to the step function as 
f{E) -^ Q{~E), so that Eqs. (7)-(9) reproduce Eqs.(6) 
and (7) in Ref. [11]. 

In our numerical calculations, we self-consistently de- 
termine Ai, (fiirj), and ^^r, for a given parameter set 
(N^, A^l, [/, 14o, Vyo, Vq, ^)- We energeticaUy compare the 
TT-phase solution with the 0-phase solution (where the 
superfluid order parameter has the same sign bwteen the 
left and right sides of the SFS-junction). The free energy 
for fixed particle numbers N^^ and N^ is given by 

F = -T^ log[l + e-^--/^] + Y, A^-^- + Ego- (10) 
i,o- (J 

III. Temperature dependence of SFS junction 

AND TT PHASE 

Figure 2 shows the temperature dependence of supefluid 
state around the SFS-junction. At T = shown in the 
upper panels, one finds that the SFS-junction and tt phase 
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Fig. 2. Calculated superfluid state in the double-well potential Vi at finite 
temperatures, in the presence of population imbalance (A'^^ = 32 > N^ = 
28). The upper, middle, and lower panels show the results for T = 0, 
0.15Tc, and 0.25Tc, respectively (where Tc = 0.4t in this calculation). 
Panels (a) and (b) show the particle densities {fii^a-) of t-spin atoms and 
^-spin atoms, respectively. Panels (c) show the local polarization s" = 
{hif) — {ni^^). Panels (d) show the spatial variation of the superfluid 
order parameter A^. We take U/t = 4. 



are realized. That is, panel (cl) shows that some of excess 
atoms are localized around the central barrier potential 
(ix ^ 0). Because of this local magnetization around the 
barrier, the superfluid order parameter A^ changes its sign 
across the junction, as shown in panel (dl), which is just 
the TT-phase. We note that, although we also obtain the 0- 
phase solution, it is energetically unfavorable, as shown 
in Fig. 3. Since the 7r-phase is known as a characteristic 
phenomenon in a ferromagnetic junction, we find that the 
localized excess atoms with pseudospin-f really behave like 
magnetic spins in the present system. 
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Fig. 3. Calculated difference of free energy AF = F-^ — Fg between 
the TT-phase (Fw) and 0-phase (Fq), as a function of temperature T. 



The TT-phase can still stably exist at finite temperatures, 
as shown in Fig. 2(d2) and Fig. 3. However, we also find in 
Fig. 2(c2) that the localized excess atoms gradually spread 
out due to thermal excitations, leading to the weakening 
of the 'magnetization' of the SFS-junction. As a result, 
the 0-phase become more stable than the vr-phase when 
T > 0.22Tc, as shown in Fig. 3. At T = 0.25Tc > 0.22rc, 
although we can still see finite magnetization of the central 
potential barrier in Fig. 2(c3), the superfluid order parame- 
ter no longer changes its sign across the junction (0-phase), 
as shown in Fig. 2(d3). 

IV. Summary 

To summarize, we have investigated the stability of tt- 
phase in a polarized superfluid Fermi gas at finite tem- 
peratures. In a model two-dimensional attractive Fermi 
Hubbard model with a double-well potential, we showed 
that the SFS-junction and vr phase can stably exist at finite 
temperatures. Our results indicate that the SFS-junction and 
TT-phase is experimentally accessible. 

In this paper, we have treated the model Fermi gas within 
the mean-field theory. In this regard, we note that pairing 
fluctuations, which are ignored in the present mean-field 
theory, are known to be crucial for ultracold Fermi gases, 
especially near the superfluid phase transition temperature 
Tc- However, since the 7r-phase is stable only far below 
Tc (T < 0.22Tc < Tc), where superfluid fluctuations 
are almost suppressed by the superfluid excitation gap, 
we expect that our results obtained in this paper is not 
drastically altered when pairing fluctuations are included. 
On the other hand, since the SFS-junction is still stable at 
T = 0.22Tc, one may need to include this strong-coupling 
effect in clarifying the temperature at which the SFS- 
junction disappears, which remains as a future problem. 
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